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Using detector tomography, we demonstrate that the detection mechanism in NbN-based su- 
perconducting single photon detectors (SSPDs) is sensitive to the total energy. In addition, we 
demonstrate a scaling law for the interchange of bias current and excitation energy. These results 
are achieved by studying multiphoton excitations in a nanodetector with a sparsity-based tomo- 
graphic method that allows factoring out of the optical absorption. We demonstrate this scaling law 
over the entire range of the experiment, from 0.9 eV to 3.8 eV. The tomographic method facilitates 
investigation of the fundamental physics of the detection events in NbN wires and demonstrates 
that nanodetectors provide a useful model system to understand more complicated wire geometries 
such as the meandering wire of a SSPD. 



I. INTRODUCTION 

Nanowire Superconducting Single Photon Detectors 
(SSPDs) [I] have high detection efficiency [2], low dark 
counts, low jitter and a broadband absorption spectrum 
[3J. This makes them suitable for many applications in- 
cluding quantum optics |1HZ], quantum key distribution 
|9j , optical time domain reflectometry [10] and inter- 
planetary communication |f f| . These detectors typically 
consist of a thin nanowire (4 nm x f 00 nm) of supercon- 
ducting material, such as NbN [J, TaN [T2], NbTiN [13], 
Nb [14] , or WSi [2], which is typically fabricated in a me- 
ander shape to cover an active area of 25-1600 /im 2 |15j . 
The absorption of a single photon in the nanowire results 
in the creation of a a region with a non-equilibrium con- 
centration of quasiparticles. When the nanowire is biased 
close to the critical current, this perturbation causes a 
transition from the superconducting to the resistive state, 
producing a voltage pulse in the external circuit. 

It has long been known that at lower bias current, 
the detector operates in a regime where multiple photons 
are necessary to break the superconductivity [T] [TD] . In 
meander-type SSPDs, this process is strongly surpressed, 
since if two photons are absorbed far apart along the 
wire, they do not constitute a single excitation and will 
not produce a detection event [T7| |T5] . 

In this publication, we use a NbN nanodetector to 
enhance these multiphoton excitations, with the aim of 
studying the physics of the detection process. We use 
a sparsity-based tomographic method, which can distin- 
guish effects at different photon numbers, to study the 
various multiphoton effects in a single experiment |19| . 
This enhancement is achieved by using a nanodetector, 
which consists of a wire with a bowtie-shaped constric- 
tion geometry (see Fig 1) [2D]. In this type of detector, 
the active area is defined by the geometric constriction, 
and all photons absorbed in this area participate in the 
detection process. 



Our results consist of three parts. 

i) We show that it is the overall energy of the excita- 
tion that determines the behaviour of the detector. An 
SSPD is therefore more properly thought of as an en- 
ergy detector rather than a single-photon detector. We 
show that the response of the device is independent of the 
wavelength of the photons that constitute the excitation 
but only dependent on the total energy of the excita- 
tion. This quantitative comparison of different experi- 
mental settings is made possible by the use of a sparsity- 
based tomographic protocol. One further advantage of 
this method is that the coupling efficiency of the photons 
to the NbN film factors out. 

ii) We demonstrate that over the range of our experi- 
ment, there is linear scaling behaviour between bias cur- 
rent and overall excitation energy. This scaling behaviour 
allows for an interchange between bias current and en- 
ergy while keeping the detection probability constant. 
All three current microscopic models for detection events 
in SSPDs predict this linear behaviour to good approxi- 
mation. 

iii) We use this scaling law to rescale all our data. We 
demonstrate that after this rescaling, all of the response 
curves superimpose, creating a photoresponse curve as a 
function of the rescaled bias current that is universal over 
the energy range of the experiment. This demonstrates 
that the detection probability in our NbN nanodetector 
is only dependent on a specific combination of bias cur- 
rent and energy. This behaviour can only be identified 
because the tomographic method separates the contribu- 
tions from different photon numbers and factors out the 
incoupling efficiency. This last result puts a constraint on 
any ab-initio theory describing the detection probability 
as a function of bias current and input energies: it must 
only depend on this given linear combination. 
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II. DETECTION EVENTS IN SSPDS 

While big strides have been taken |21lf2"4"] in under- 
standing the fundamental physics of these detectors, 
many details of the detection mechanism in such detec- 
tors are still unknown. After the absorption of a pho- 
ton, it is thought that the resulting high-energy elec- 
tron destroys Cooper pairs that carry the bias current, 
producing a cloud of quasi-particles [SB]. In some 
models, this breakdown is assisted by a thermal fluctu- 
ation process, such as a vortex crossing, which adds an 
extra perturbtion to the system. This process results 
in a breakdown of the superconductivity, resulting in a 
normal cross-section. After such a resistive barrier has 
formed, the kinetic inductance of the device drives Joule 
heating in the normal state area [57], causing it to grow. 
After that, the current drops to a negligible level [2"8] 
and the interplay between the cooling of the device and 
the restoration of the current determines whether the de- 
vice returns to its previous state, ready to detect another 
photon. 

Currently, there is no consensus on a microscopic 
model for the detection event in superconducting single 
photon detectors. Below, we briefly introduce three mi- 
croscopic models for detection events in the SSPD from 
literature. We refer to these three as the normal-core 
hotspot model, the diffusion hotspot model and the fluc- 
tuation model. Each of these models predicts that the en- 
ergy and current are exchangeable through a given scal- 
ing law, and their prediction will be compared with the 
experimental data below. 

The normal-core hotspot model was introduced in the 
original papers reporting photodetection with SSPDs 
[TJ [29] . In this model, it is assumed that the photon ab- 
sorption creates a normal core inside the material. Cur- 
rent is then diverted around this core. If the current 
locally exceeds the critical current, superconductivity is 
destroyed and a normal state slab is created, resulting 
in a detection event. In this model, the current required 
to pinch off the entire channel for a given input energy 
goes as the square root of the energy, since the hotspot is 
assumed to be a cylindrical object inside the wire. For a 
given energy, the bias current It, needed to achieve max- 
imal detection efficiency is then given by: 

E=^{\-hlh)\ (1) 

where E is the energy of the photon, is the bias current 
and I c is the critical current, w is the width of the wire 
and C is a scaling constant, which is defined in this way 
for consistency with previous work [30J. 

In the diffusion-hotspot model, which was introduced 
later as a refinement of the original normal-core hotspot 
model, the role of the diffusion of quasiparticles is taken 
into account, as well as the reduction of the critical cur- 
rent due to the quasiparticles |2T| [25] . In this model, the 
relevant lengthscale is given by the diffusion length over a 
time characteristic for the cascade of quasiparticles. This 



expression, which was first derived in Ref. [25], predicts 
linear scaling between bias current and cutoff energy: 

E = E (l-I b /I c ), (2) 

where Eq is an energy scale [26 . 

Fluctuation model The previous two models predict 
a sharp cutoff of the detection efficiency as a function 
of photon energy, which is not observed in experiments. 
In order to explain the observed detection probability 
beyond the cutoff energy, fluctuation-assisted detection 
models have been proposed [22--25J. All of these have in 
common that the role of the photon in the detection pro- 
cess is to depress the superconducting gap. Subsequently, 
a thermally activated fluctuation occurs, such as the de- 
pairing of a vortex-antivortex pair (VAP) [31 in the su- 
perconductor or the crossing of a single vortex. This fluc- 
tuation must overcome an energy barrier E(A,Ib) |24| . 
Expressions for such energy barriers typically contain a 
gap-dependent energy scale and a current-dependent ge- 
ometric factor [26 . The specifics of the geometric fac- 
tor depends on the precise fluctuation process. Both for 
the VAP model and for the single- vortex crossing model, 
we can linearize the current-dependence of the geometric 
factor over the range of currents that was used in the 
experimernt. We obtain: 

A=(A-aVE)(I -f3h), (3) 

where the constants Iq and j3 are known from the lin- 
earization of the geometric factor, and A and a are the 
experimental fit parameters. Such a model predicts a 
hyperbolic interchange between energy and current. 

III. EXPERIMENT 

All experiments in this work were performed on a nan- 
odetector (see Fig 1). The nanodetector consists of 4 nm 
thin NbN wire on a GaAs substrate, shaped into a 150 
nm wide bowtie geometry. The device was fabricated via 
a combination of DC magnetron sputtering [32], electron- 
beam litography, reactive ion etching and evaporation of 
the metal contacts [20]. In previous work [191 120] . it 
was shown that such a detector has multiphoton regimes 
based on the bias current. The physical mechanism be- 
hind these multiphoton regimes is that at relatively low 
bias currents, multiple photons are required to supply 
a sufficient perturbation for the superconductivity to be 
broken. 

The device was cooled in a two-stage pulse-tube / 
Joule-Thompson cryocooler to a temperature of approx- 
imately 1.2K. The nanodetector was illuminated using a 
lensed fiber mounted on cryogenic nanomanipulators. At 
this temperature, the overall system detection efficiency 
for single photons was 1.5 x 10~ 4 around our working 
point at lb — 20 //A (I c = 29 fj,A) . This low efficiency is 
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Figure 1: False-color SEM image of the detector. The active 
part of the detector is the narrow bridge in the centre of the 
image. The blue parts represent the thin layer of NbN, the 
red parts are the GaAs substrate. The scale bar has a length 
of 1 fim. 

attributable to the mismatch between the device active 
area and the size of the illumination spot. The device 
was operated in a voltage bias regime, using a Yokogawa 
GS200 in series with a 10 resistor. The detector was 
biased through the DC port of a bias tee, and the RF 
clicks were amplified in a 45 dB amplifier chain. 

The device was illuminated with a Fianium Supercon- 
tinuum laser, whose pulse duration was specified to be 7 
ps. It is crucial for this experiment that the pulse dura- 
tion is shorter than the lifetime of an excitation, which 
is estimated to be several tens of picoseconds [33 35 ; for 
longer pulse durations, it is possible to have a pulse which 
produces two excitations which are far enough apart in 
time as to not produce a single multiphoton excitation. 
The device was illuminated with narrowband light at 
wavelengths of 1000 nm, 1300 nm and 1500 nm (AA= 
10 nm). In our experiment, we vary the intensity and 
wavelength of the input light, at various bias currents. 
At each of these settings, we record the count rates in a 
0.1 s time window and repeat the experiment 10 times 
per measurement setting. 

IV. TOMOGRAPHY OF MULTIPHOTON 
EXCITATIONS 
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Figure 2: Illustration of the tomographic protocol. The black 
squares indicate the measured count rate as a function of in- 
put power, at A = 1500 nm and h = 17 //A. The red and 
green lines show the contribution to the count rate of single 
photons and photon pairs, respectively. The blue line shows 
the contribution of higher numbers of photons. The black 
line shows the sum of all the photon contributions, indicating 
that our tomographic reconstruction succesfully reproduces 
the observed count rates. From this fit, we reconstruct the 
set of detection probabilities p n and r\ which describes the 
behaviour of the detector. 

This model was introduced in a previous publication [TU] , 
and we give here a brief summary of the main features. 
In this protocol, it is assumed that the number of param- 
eters that describes the detector is small, and that most 
of the behaviour of the count rates is dominated by low 
incoupling efficiency. Using model selection rules it 
is possible to find the description of the detector that 
requires the least number of parameters to be fully con- 
sistent with the experimental data , thus invoking the 
sparsity assumption. In particular, we model the click 
probability R c u ck by: 



In order to distinguish the effects of the various pho- 
ton numbers in the laser pulses, we make use of a to- 
mographic protocol. We illuminate the detector with a 
range of coherent states, and record the click probabil- 
ity Rdick- Each illumination intensity probes the detec- 
tor with a different linear combination of photon num- 
bers, introducing different combinations of multiphoton 
excitations. Since we have many more experiments than 
photon numbers that are of interest, we obtain an over- 
complete set of equations, which can be inverted using 
suitable techniques for such problems |36H38| . 

Since the physical size of the detector is much smaller 
than the wavelength, the incoupling must necessarily be 
very inefficient due to the diffraction-limited size of the 
illumination spot. To perform tomography in this geom- 
etry, we have established a sparse tomographic protocol. 



Rcuc = 1 - e-i» V(l - Pn )^-, (4) 

^— ' TV. 

Tl = 

where rj is the incoupling efficiency, and N is the mean 
photon number of the incident coherent state. The mean- 
ing of the set of p n is as follows: the pi are the param- 
eters that describe all nonlinear (multiphoton) effects in 
the detector. They are therefore the quantity of interest 
in the further analysis of the behaviour of the detector. 
The linear efficiency rj quantifies the linear effects in the 
detector, i.e. the probability that each individual pho- 
ton participates in the detection process. In making this 
separation, we have used the fact that an attenuated co- 
herent state is still a coherent state, with reduced mean 
photon number. This is reflected by the fact that in Eq. 
4, only the combination r/N occurs. 
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We attribute the linear loss quantified by r) to the com- 
bination of coupling of the beam to the detector area and 
absorption in the wires • Each p n indicates the prob- 
ability that the detector clicks, given that n photons are 
absorbed in the active area and participate in the detec- 
tion event. The probability of this occuring is assumed 
to be independent for each photon, and is given by the 
linear efficiency ry. The separation between 77 and {p n } 
is therefore one between linear and nonlinear effects, re- 
spectively. 

The infinite sum in Eq. 4 is truncated at some n max 
by our model selection rule: we select the model with 
minimal reduced x 2 . This truncation enables us to re- 
cover both 77 and the set of nonlinear parameters p n from 
a single experiment. 

Fig. 2 illustrates this protocol as applied to a single 
experimental run for a given bias current. We vary the 
incident power, observe the click probability, and apply 
the procedure outlined in this section to find the contri- 
butions from the various multiphoton excitations. The 
black squares indicate the measured count probability, 
approaching 1 as the detector saturates. The red (solid), 
green (dashed) and blue (dotted) lines indicate the contri- 
bution from one photon, two photons and higher photon 
numbers, respectively. The remainder term, which con- 
tains the limit of high photon numbers, is not used in the 
results presented below. Fig. 2 shows that there is a low- 
power regime where the detector responds primarily to 
single-photon excitations, an intermediate regime where 
the detector responds to two-photon excitations, and a 
third regime where the detector is approaching satura- 
tion. The fact that the photon regimes split in this way 
with input power enables us to recover them seperately. 
Furthermore, since the linear efficiency 77 only rescales 
the effective incident photon number but does not alter 
its shape (corresponding to a simple shift in Fig 2), we 
are also able to distinguish finite incoupling effects from 
effects due to multiphoton excitations. 



V. RESULTS 

Figures. 3-5 present our main results. Fig 3 shows 
the reconstructed detection probabilities p n , as a func- 
tion of bias current and for different wavelengths. For 
each wavelength and current, we independently perform 
the tomographic procedure outlined in section IV and 
obtain a full set of nonlinear parameters p n .We observe 
that as the current is lowered, the detector makes a tran- 
sition from being a one-photon detector to a two-photon 
detector, and so on. Furthermore, we observe that the 
response curves at different photon numbers and wave- 
lengths have the same shape. We note that as the ex- 
citation energy becomes higher and the photon number 
larger, the points on our curves become more scattered, 
indicating that the tomography procedure becomes less 
accurate. 

Fig. 4 demonstrates that there is a scaling law be- 
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Figure 3: Current dependence of the nonlinear parameters 
p n , as a function of wavelength and photon number. The 
probability p n of the detector clicking at a given wavelength 
and photon number is plotted as a function of the current. 
The plots are color-coded by wavelength. The shape of the 
symbols indicates the photon number (see legend). The con- 
necting lines are a guide to the eye. The dotted line indicates 
the threshold level used to obtain Fig. 4. 
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Figure 4: Scaling law for the nanodectector. From the data 
in Fig 3, we find all points that have p n (E,Ib) = 0.1 (indi- 
cated by the dotted line in that figure), where E is the overall 
excitation energy. In this graph, we plot the values of It and 
E that satisfy this condition. This graph shows that bias cur- 
rent and overall excitation energy have a linear dependence. 
Furthermore, the fact that points at various photon numbers 
all fall on the same line demonstrates that the nanodetector 
is only sensitive to the overall energy of the excitation. The 
three lines show the fits of the three microscopic models to 
the data. 
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VI. DISCUSSION 
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Figure 5: Universal response curve for the nanodetector. To 
obtain these curves, we rescale the curves reported in Fig. 3 
by the scaling law demonstrated in Fig. 4. 



tween bias current and overall excitation energy. In or- 
der to obtain this figure, we took a surface of constant 
p n (E,Ib) = 0.1 in Fig. 3 (indicated by a dotted line), 
and plotted the bias current at which the detector has 
10% probability of responding to an energy E, where 
E = nhv is the total energy of the n photons absorbed 
by the detector. The three curves in Fig. 4 show the best 
fit of the three microscopic models to the data. Over the 
range of the experiment, the predictions for the three 
models are to good approximation equivalent: they all 
predict the observed linear interchange between photon 
energy and bias current. Only below E < 0.5 eV do the 
three models differ significantly. Fig. 4 shows that there 
is a scaling law between current and energy: in our de- 
tector, they are interchangebale at a rate 7 of 7 = -2.9 
± 0.1 /j,A/eV (= -1.8 x 10 13 Wb-Hn SI units). 

In Fig 5, we apply the scaling law and find that all the 
curves of detection probability as a function of rescaled 
bias current superimpose over 5 orders of magnitude in 
the detection probability. This shows that the photore- 
sponse of our detector depends only on this specific com- 
bination of bias current and excitation energy. We stress 
that this universal curve can only be obtained through 
detector tomography that allows separation of the effects 
of multiphoton excitation and finite linear efficiency. 

Since the curve is universal over 5 orders of magnitude 
in the detection probability, we have demonstrated that 
our results are independent of the arbitrary choice of the 
10% criterion. The criterion only matters for the ac- 
curacy with which the curves can be superimposed: we 
find from theoretical simulations that the tomographic 
reconstruction is most accurate between p n — 0.1 and 
p n = 10 -5 . This justifies the choice of our criterion. 



In this section, we first compare our experimental 
method with that of previous studies on SSPDs. Then, 
we discuss our three experimental findings: energy de- 
tection, linear interchange between bias and energy and 
the universal detection curve. Lastly, we discuss the phe- 
nomenon of dark counts in our detector. 

Comparison with previous work Our present method 
has three advantages over methods used in previous stud- 
ies, in which SSPDs were studied by varying the wave- 
length of the incident light [25, 26 . First of all, through 
a tomographic method, we obtain a full quantitative de- 
scription of our detector. This description allows com- 
parison of different detection settings, enabling us to find 
the scaling law reported in Section IV. Secondly, our to- 
mographic protocol can distinguish various multiphoton 
excitations, which enables us to extend the results from 
previous work to the multiphoton regime. Thirdly, since 
the linear absorptance is factored out by the fitting pro- 
cedure, the derived detector response parameters p n are 
independent of the coupling efficiency and absorptance. 
In particular, the wavelength dependence of the optics 
inside the cryostat, of the absorption into the NbN and 
of any associated cavity structure [ID] is removed. This 
strongly facilitates quantitative comparison between dif- 
ferent experimental settings. 

Multiphoton excitation has the very practical advan- 
tage that the bandwidth of energy excitations which is 
offered can be extended by a factor equal to the num- 
ber of photons in the highest excitation (in our case, 4). 
This has applications in the situation where light of a 
particular wavelength is difficult to couple onto a cryo- 
genic sample. In this way, we extend the available en- 
ergy range of excitations for any given cryogenic setup. 
In particular, the present work opens up the possibility 
of studying NbN detector behaviour in an energy range 
which corresponds to the near and medium UV range. 

The tomographic method does not require a model of 
the device, which makes it useful as a tool for investigat- 
ing the working principle of any type of detector. SSPDs 
are an important case study for the application of the 
tomographic method, since they possess both nontrivial 
detection statistics |17[fT8]. and their working mechanism 
is not fully understood yet. 

Energy detection. The fact that only the overall energy 
of the excitation determines the response of the photode- 
tector can be interpreted in terms of the cascade process 
that generates the initial excitation. This process, which 
is thought to involve both electrons and phonons in the 
film, and in which the mutual exchange of energy between 
the electron and phonon subsystem plays a key role, is 
still poorly understood. In the present work, we probe 
this cascade process with different initial excitations, and 
show that it is only the overall energy which determines 
the total number of quasiparticles which are produced 
at the superconducting band-edge. The fact that four 
excitations of a quarter of the energy produce the same 
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number of QP as a single excitation with the full energy 
is evidence of the fact that the conversion efficiency by 
which the energy of the first QP is distributed over many 
others is independent of the initial energy. 

In the normal-core hotspot model, multiphoton events 
can only occur if the various normal cores are created at 
the cross-section of the wire. Given the relative size of our 
detector and the size of the normal cores in this hotspot 
model, we would expect to see reduced count rates at 
high photon numbers due to the low probability of several 
photons being absorbed at exactly the same point in our 
detector. The fact that we do not observe such reduced 
count rates is evidence the in our experiment, the size of 
the hotspot is comparable to the size of our detector. 

Linear interchange We apply the three models, and 
compare the results with the values from literature. 
For the normal-core hotspot model, we find C — 46 
eV- 1 / 2 / nm, which should be compared to the values 
of C = 11 - 20 eV" 1 / 2 / nm found in other experiments 
[3]- For the diffusion hotspot model, we apply the ex- 
pression from Ref . [26_ , to find a theoretical value of 7 
= -2.5 /iA/eV for our sample and 7 = -3.5 fiA/eV for 
the samples in that reference, which should be compared 
with the value of 7 = -2.9 /iA/eV obtained experimen- 
tally. We note, however, that this result depends linearly 
on the assumed conversion efficiency of the initial ex- 
citation into quasiparticles, which varies strongly from 
detector to detector. For the fluctuation model, we find 
a = 2.2 x 10 -4 \/eV, which should be compared to a lit- 
erature value of a = 6 x 10~ 4 \/cV, for the experiment 
reported in Ref. |31| . 

We note that in all three models, there is one overall 
scale parameter which is meaningless, because its value 
depends on the precise choice of the threshold criterion 
in Fig 3. In the two hotspot models, this introduces 
an offset current. In the fluctuation model, the energy 
scaling constant A depends on the choice of the threshold 
criterion. 

Universal detection curve By demonstrating a univer- 
sal detection curve, we have shown that the various mod- 
els for photodetection do not just predict the scaling be- 
haviour for the plateau region of the detection efficiency, 
but for the entire detection curve. Such quantitative 
comparison would not be possible without the use of the 
tomographic method. 

In a previous publication |35| , we have introduced the 
notion of the nonlinear response function (NRF) t](Ii,, C), 
which measures the instantaneous click probability, given 
that a bias current of 2b is present, and that there are C 
quasiparticles in the detector. The overall click proba- 
bility is then given by R = J r)(t)I(t)dt, where I is the 
instantaneous intensity. This function may be probed 
through various means, such as measuring the autocor- 
relation in a pump-probe experiment. The description in 
terms of a NRF is well-matched to a tomographic exper- 
iment, as both are model-independent descriptions. 

The holy grail of tomographic research on SSPDs 
would be to find the instantaneous click probability as a 



function of the number of quasiparticles present at that 
instant. In the present experiment, we have achieved a 
step towards this goal: we have demonstrated the NRF 
to be of the form 77 (If, + ^E) over the energy range of the 
experiment. 

Dark counts We now turn to the phenomenon of dark 
counts. Extrapolating the linear scaling law from Fig. 4 
to E = yields a current of 19 /iA. However, the actual 
dark count rate at this current was negligibly small; we 
only observe appreciable dark counts around the critical 
current of 29 /iA. The same issue holds for the extrapola- 
tions from the other two models. Therefore, our scaling 
law does not hold for dark counts. For this, we propose 
the following explanation: due to the geometry of our 
detector, there will always be a vortex pair present at 
the corners of the detector. When a photon is absorbed, 
the vortex pair is depinned, and the energy contained 
in it is added to the excitation. However, when there is 
no photon present, the vortex pair must tunnel through 
an energy barrier, producing a much lower dark count 
rate. The energy barrier for a single vortex to cross our 
detector is of the order of 20 meV in our sample [26 . 

The anomalous behaviour of dark counts is a reminder 
of the danger of assuming a detection model, further 
demonstrating the relevance of our tomographic method. 
In particular in this case, the tomographic method may 
give the first hints of substantial differences in detection 
mechanism between dark counts and light counts. While 
we have proposed an explanation for this behaviour in 
terms of vortex pairs which is consistent with the data, 
we note that the nature of dark counts is still open to 
debate 



VII. CONCLUSIONS 



We have studied the physics of photodetection in a su- 
perconducting single photon detector. We have shown 
three main results. First, we have shown that the detec- 
tion is based on the overall energy of the excitation. Sec- 
ond, that there is a scaling law between overall excitation 
energy and bias current. Third, that there is a universal 
response curve that depends only on a given combination 
of bias current and excitation energy. Thereby, we have 
shown that the known behaviour of the detector extends 
into the multiphoton range. These results demonstrate 
that the tomographic method is a useful tool for inves- 
tigating the fundamenal physics of detection events in 
NbN SSPDs. Our technique extends the available range 
of energy excitations beyond that of a given cryogenic 
setup, and we have shown that the nanodetector is a 
useful model of a single cross-section of a meander-type 
SSPD. 
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